How many young star clusters exist in the Galactic center? 
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ABSTRACT 

We study the evolution and observability of young compact star clusters within 
~ 200 pc of the Galactic center. Calculations are performed using direct TV-body 
integration on the GRAPE-4, including the effects of both stellar and binary 
evolution and the external influence of the Galaxy. The results of these detailed 
calculations are used to calibrate a simplified model applicable over a wider 
range of cluster initial conditions. We find that clusters within 200 pc of the 
Galactic center dissolve within ~ 70 Myr. However, their projected densities 
drop below the background density in the direction of the Galactic center 
within ~ 20 Myr, effectively making these clusters undetectable after that time. 
Clusters farther from the Galactic center but at the same projected distance are 
more strongly affected by this selection effect, and may go undetected for their 
entire lifetimes. Based on these findings, wc conclude that the region within 200 
pc of the Galactic center could easily harbor some 50 clusters with properties 
similar to those of the Arches or the Quintuplet systems. 

1. Introduction 

Two young compact star clusters have been observed within a few tens of parsecs 
of the Galactic center: the Arches cluster (Object 17, Nagata ct al. 1995) and the 
Quintuplet cluster (AFGL2004, Nagata ct al. 1990; Okuda et al. 1990), for which excellent 
observational data are available. In terms of structural parameters — size, mass, density 
profile — and ages, these systems may represent the Galactic counterparts to NGC2070 
(R136), the central star cluster in the 30Doradus region in the Large Magellanic Could 
(Massey & Hunter 1998). The Arches and Quintuplet clusters he behind thick layers of 
absorbing material, hinting that many more such systems may exist. Recently, Dutra & 
Bica (2000) have reported from the 2MASS survey a total of 58 star cluster candidates 
within ~ 600 pc (in projection) of the Galactic center. 

A number of important questions make these clusters worthy of detailed study, among 
them: (1) How are such clusters related to globular clusters? (2) How do they contribute 
to the total star formation rate in the Galaxy? (3) Are their mass functions in reality 
intrinsically fiat, as is suggested by observations? (4) How far are these clusters from the 
Galactic center? (5) How many are hidden, still waiting to be discovered? In this Letter we 
summarize the results of a series of A-body simulations of young compact star clusters in 
the vicinity of the Galactic center, and address specifically the last item on this list. 
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We found that the hfetimes of our model clusters depend very sensitively on their 
distances from the Galactic center. This is mainly due to the larger size of tidally limited 
clusters lying farther from the Galactic center, resulting in longer relaxation times and 
therefore longer lifetimes. The majority of our models are visible only for the first part of 
their lifetimes, and are likely to be indistinguishable from the stellar background at later 
times. We find that the true number of young compact star clusters within 200 pc of the 
Galactic center is at least 10 but could easily exceed 50. 

A more comprehensive paper, exploring all of these questions in more detail and 
presenting the results of an extensive parameter-space study, is in preparation (Portegies 
Zwart et al., 2000b). 

2. Initial conditions 

We study the evolution of our model clusters by following the equations of motion of 
all stars by direct A'"-body integration, at the same time taking into account the internal 
evolution of both stars and binary systems. The "Starlab" software environment within 
which this work was performed is described in detail by Portegies Zwart et al. (2000a; 
see also pittp : //www, sns . ias . edu/^starlab|) . The special-purpose GRAPE-4 system 
(Makino et al. 1997) was used to accelerate the computation of gravitational forces between 
stars. 

Observed parameters for the Arches and Quintuplet clusters are listed in Table |I]. 
These clusters have masses of ~ IO'^Mq and are extremely compact, with half-mass radii 
'"hm ~ 1 pc (Figer, McLean & Morris 1999). The projected distance from the Arches to the 
Galactic center is about 34 pc; the Quintuplet cluster lies somewhat farther out, at ~50pc. 

Our calculations start with 12k (12288) stars at zero age. We choose stellar masses 
between 0.1 M© and 100 M©, distributed according to the mass function suggested for the 
Solar neighborhood by Scalo (1986). The median mass of this mass function is about 0.3 M©; 
the mean mass is (m) ^ 0.6 M©. The initial mass of each model is therefore ~ 7500 M©. 
Initially all stars are single, although some binaries do form dynamically via three-body 
encounters, in which one star carries away sufficient energy and angular momentum to allow 
two others to become bound. We adopt three standard distances from the Galactic center, 
34 pc, 90 pc and 150 pc. The initial density profiles and velocity dispersions for our models 
are generated from anisotropic Heggie & Ramamani (1995) models with Wq = 4. At birth, 
the clusters are assumed to precisely fill their Jacobi surfaces ("Roche lobes") in the tidal 
field of the Galaxy, and are taken to move in circular orbits around the Galactic center. For 
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Table 1: Observed parameters for the Arches and the Quintuplet clusters. Columns list 
cluster name, reference, age, mass, projected distance to the Galactic center, tidal radius 
('^tide), and half mass radius (rhm)- The final column presents an estimate of the density 
within the half mass radius. 



Name 


ref 


Age 


M 




?"tide 




Phm 






[Myr] 


[103 Mq] 




-[pc] 




[IO^Mq/pc^] 


Arches 


a 


1-2 


12-50 


30 


1 


0.2 


0.6 - 2.6 


Quintuplet 


b 


3-5 


10-16 


50 


1 


~ 0.5 


0.08 - 0.13 



References: a) Brandl et al. (1996); Campbell et al. (1992); Massey & Hunter (1998). b) 
Figer et al. (1999); 



a circular orbit in the plane of the Galaxy the distance from the center of the star cluster 
to the first Lagrangian point (the Jacobi radius) is approximated by 

Here M is the mass of the star cluster and the factor two is a correction factor which 
depends on the the density profile; strictly speaking the factor of two is correct only in the 
case Mcai oc r. (The Jacobi radius is computed consistently with the adopted tidal field in 
our simulations.) Table ^ presents an overview of our model initial conditions. 

Table 2: Overview of initial conditions for our model calculations. Each calculation is 
performed three times. From left to right the columns list the model name, the distance 
to the Galactic center, the initial King parameter Wo, the initial tidal- and half mass 
relaxation times, half-mass crossing time, core radius, half-mass radius and distance to the 
first Lagrangian point in the tidal field of the Galaxy, the time of core collapse, and the time 
at which the cluster mass drops below 1% of its initial value. 



Model 




Wo 


^rxt ^rxh 




'"core 










[pc] 




— [Myr]— 


[kyr] 




- [pc] - 




— [Myr]— 


R34W4 


34 


4 


53 3.2 


27 


0.05 


0.117 


0.77 


0.8 12.7 


R90W4 


90 


4 


134 8.1 


68 


0.09 


0.218 


1.42 


1.2 32.6 


R150W4 


150 


4 


218 13 


110 


0.14 


0.301 


1.97 


2.0 53.4 



The mass of the Galaxy within the clusters' orbit at a distance tqc ( ~ 200 pc) is taken 
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to be (Mezger et al. 1999) 

MGai(rGc) = 4.25 X 10^ (^^^ [M©]. (2) 

This mass distribution determines the strength and geometry of the local Galactic tidal 
field (for details see Portegies Zwart et al. 2000b). The evolution of the cluster is followed 
using the Starlab kira TV-body integrator and the SeBa binary evolution program (Portegies 
Zwart et al. 2000b). For each selected distance to the Galactic center we carried out three 
calculations. One series of runs was carried out with identical initial realizations of the 
A^-body system (stellar masses, positions and velocities, with a total mass of 7432 M©). The 
same initial model can be used at several galactocentric distances because the shape of the 
zero- velocity surface does not depend sensitively on distance to the Galactic center. Table ^ 
gives results of these calculations. For each galactocentric distance we also performed two 
additional calculations (for a total of 9 runs) with different initial realizations of the N-hodj 
systems. These calculations were performed to study the uncertainties in cluster lifetimes, 
and to ascertain the reproducilility of our results. The calculations with different initial 
realizations produced roughly 10% spreads in core collapse times (tec) and cluster lifetimes 
(tend)- For reasons of economy, stars were removed from all iV-body calculations when they 
exceeded a distance of 3 tli from the cluster center. 



3. Results 

Figure |1| shows the evolution of cluster mass and number of stars for the models 
listed in Table H which began with identical initial conditions but different galactocentric 
distances. Perhaps not surprisingly, clusters located farther from the Galactic center live 
considerably longer than those closer in. The longer lifetime of the more distant cluster 
is mainly a consequence of its longer relaxation time. Scaling the relaxation time at the 
tidal radius of model R34W4 to a distance of 150 pc results in a lifetime of ~ 52.2 Myr 
(= 12.7Myr x 218/53), which is slightly smaller than the ~ 53.4 Myr lifetime of the models 
R150W4. Mass loss from stellar evolution, which is more prominent in model R150W4 
because the cluster survives longer, seems to be a minor factor in driving the evolution. 

The number of stars in each model (dashed lines, see Figure 1) decreases more rapidly 
than the total mass (solid lines). Thus the mean mass of the stars within the cluster 
increases gradually with time. 
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Fig. 1. — Evolution of the total mass M and number of stars (solid and dashed lines 
respectively) within the critical zero- velocity surface of selected models at rcc = 34 pc (left), 
'"GC = 90 pc (middle) and rcc = 150 pc (right). The results of the runs with identical initial 
realizations are presented. All quantities are normalized to their initial values. 
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4. Discussion 

Although clusters like the Arches and Quintuplet systems are very compact, it may 
still be hard to see them near the Galactic center because the projected foreground and 
background stellar density is so high. Integration of the local stellar density (obtained by 
differentiating Eq. ^ along the line of sight then gives the surface density. Portegies Zwart 
et al. (2000b) perform this calculation numerically and arrive at a surface density at 34 pc 
of about 3000 MqPc"^. While the contrast in the surface density of the cluster relative to 
that of the background is perhaps an oversimplified measure of the cluster's observability, 
the results of this simple comparison are quite instructive, and a more comprehensive 
consideration of luminosity density leads to essentially similar overall conclusions. 




t [Myr] 



Fig. 2. — Evolution of the surface density within the projected half-mass radius for the 
models at rcc = 34 pc (dotted line), at rcc = 90 pc (dashed line) and at rcc = 150 pc (solid 
line). The horizontal dotted line gives the integrated background density at a projected 
distance of 34 pc from the Galactic center. The results of the runs with identical initial 
realizations are presented. The two error bars give the observed surface densities of the 
Arches (left) and the Quintuplet (right) clusters. 

Figure|^ shows the evolution of the surface density within the projected half mass radius 
for models R34W4 (dots), R90W4 (dashes) and R150W4 (solid). The two points with error 
bars indicate the projected half-mass densities for the Arches (left) and Quintuplet (right) 
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clusters. The horizontal dotted line gives the background surface density at a projected 
distance of 34 pc from the Galactic center. The projected densities of the two observed 
clusters are between 3 (for the Quintuplet) and 50 (for Arches) times higher than the 
background; clusters with densities below a few times the projected background stellar 
density may well remain unnoticed. 



4.1. A simple model 

A simplified model for the evolution of these star clusters may be constructed as 
follows. The initial relaxation time at the tidal radius may be calculated using Eq. ^ and 
Spitzer's (1987) expression as 

/ \ 0.9 

t„t^2.19[^j [Myr]. (3) 

The constant is obtained by substitution of the appropriate units. 

The mass of the cluster decreases almost linearly in time (see Figure |l]) as 

M = Mo(l-yy (4) 

Here r = t/trxt and Tc ^ 0.29trxt is the age at which the cluster dissolves in the tidal field of 
the Galaxy. The projected surface density within the half-mass radius of the cluster is 

Mhm /■r\ 
E = = (5) 

Here, Mhm — 0.65M is the mass contained within the projected half-mass radius, and 
Wo = Thm/''"Li depends on the density profile, but is always smaller than unity. For a King 
model with Wq = 4, we find Wo — 0.16. Substitution of Eqs.|I] and into Eq.|^ gives 

So ^ 7.0 X 10^ (^j^ j (6) 

The projected surface density E decreases with time because the cluster mass decreases 
and the half-mass radius of the cluster increases. Substitution of Eq. ^ into Eq. ^ gives: 

Eo V TcJ \ Thm / 

The half-mass radius increases by about a factor of two during the first half mass 
relaxation time and remains roughly constant thereafter. In our simple model we implement 
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this by allowing ~ '''hm.o to increase by a factor two in the first half mass relaxation 
time and to remain constant at later times. The fact that the half-mass radius remains 
roughly constant at late times, and does not decrease as M^/^ as would be expected for a 
tidally limited system, is a consequence of our use of the total TV-body mass, rather than 
the mass within the Jacobi surface (see Figure P, in determining both rhm and S. The 
resulting surface density evolution agrees very well with our A^-body calculations. 




Fig. 3. — Projected cluster density as a function of distance from the Galactic center and 
time. Gray shading indicates projected surface density; darker shades indicate higher density. 
Solid lines indicate where the surface density of the cluster equals 1, 2, 4, and 10 times the 
background density. The background surface density is computed for 34 pc. The two holes 
to the lower left indicate the locations of the Arches (left) and Quintuplet (right) clusters 
on the figure. 

Figure]^ shows the evolution of the projected density as a function of distance from 
the Galactic center for our simple model. Clusters close to the Galactic center are compact 
enough to be easily visible (dark shades) for a large fraction of their lifetimes, but they 
dissolve quickly. Clusters farther from the Galactic center live much longer, but their surface 
densities are generally low, lying well below the background for most of their lifetimes. The 
two well-observed clusters both lie in the lower left corner (high density region) of Figure^ 
The Arches and Quintuplet clusters inhabit only a small portion of the available parameter 
space. However, the region they populate is most favorable for finding clusters, because the 
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projected surface densities of such clusters are high. 

We estimate the number of clusters like the Arches and the Quintuplet based on the 
fact that both clusters lie at the lower left corner (youngest and most compact) in Figure |[ 
We assume that all clusters with ages less than ~ 5Myr and within 50 pc of the Galactic 
center have been found by observers, and all outside this region are as yet undetected. 
We further assume that clusters populate the triangular area in Figure ^ more or less 
uniformly — probably not a bad assumption, since the Galactic mass is roughly proportional 
to radius (see Equation||). With these assumptions, the number of yet-to-be-found clusters 
is (ix60Myrx200pc)/(5Myrx50pc)=24 times more than the number of known clusters. 
Thus, we expect that the total number of clusters in this region would be around 50. 

A conservative lower limit to the number of hidden clusters may be obtained using the 
same technique, but adopting the observed projected density of the Quintuplet system as 
the limiting contrast at which such clusters can be discovered. The projected density of the 
Quintuplet exceeds the background density by about a factor three. Figure^ then suggests 
that about 20% of the available surface area harbors visible cluster. A lower limit to the 
total number of clusters within 200 pc of the Galactic center would then be around 10. 

These estimates are in excellent agreement with the results of Dutra & Bica (2000) 
mentioned earlier, although the above reasoning suggests that even these 58 candidate 
clusters may represent only a small fraction of the number actually present. 

Finally, we note that a population of 100 clusters with masses of 10^ M© each and a 
maximum lifetime of lO^Myr implies a star formation rate of 0.01 per year, enough to 
build up the entire bulge of 10^ Mq within the 10 Gyr age of the Galaxy. 

SPZ is grateful to the Institute for Advanced Study, Drexel University and Tokyo 
University for their hospitality and the use of their GRAPE-4 hardware. This work was 
supported by NASA through Hubble Fellowship grant HF-01112.01-98A awarded by the 
Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, by the Research for the Future Program of Japan Society for 
the Promotion of Science (JSPS-RFTP97P01102) and by NASA ATP grants NAG5-6964 
and NAG5-9264. Part of this letter was written while SPZ, SM and PH were visiting 
the American Museum of Natural History. They acknowledge the hospitality of their 
astrophysics department and visualization group. 



- 11 - 



REFERENCES 

Brandl, B., Sams, B. J., Bertoldi, F., Eckart, A., Genzel, R., Drapatz, S., Hofmann, R., 
Loewe, M., Quirrenbach, A. 1996, ApJ, 466, 254 

Campbell, B., Hunter, D. A., Holtzman, J. A., Lauer, T. R., Shayer, E. J., Code, A., Faber, 
S. M., Groth, E. J., Light, R. M., Lynds, R., O'Neil, E. J., J., Westphal, J. A. 1992, 
AJ, 104, 1721 

Dutra, CM., & Bica, E., 2000, A&A, 359, L9 

Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. M., McLean, L S. 1999, ApJ, 
525, 750 

Figer, D. F., McLean, I. S., Morris, M. 1999, ApJ, 514, 202 
Heggie, D. C, Ramamani, N. 1995, MNRAS, 272, 317 
Makino, J., Taiji, M., Ebisuzaki, T., Sugimoto, D. 1997, ApJ, 480, 432 
Massey, R, Hunter, D. A. 1998, ApJ, 493, 180 

Mezger, R C, Zylka, R., Philipp, S., Launhardt, R. 1999, A&A, 348, 457 

Nagata, T., Woodward, C. E., Shure, M., Pipher, J. L., Okuda, H. 1990, ApJ, 351, 83 

Nagata, T., Woodward, C. E., Shure, M., Kobayashi, N. 1995, AJ, 109, 1676 

Okuda, H., Shibai, H., Nakagawa, T., Matsuhara, H., Kobayashi, Y., Kaifu, N., Nagata, T., 
Gatley, L, Geballe, T. R. 1990, ApJ, 351, 89 

Portegies Zwart, S. F., Verbunt, F. 1996, A&A, 309, 179 

Portegies Zwart, S. F., McMillan, F., Makino, J., Hut, P. 2000a, MNRAS in press 

(|stro^pE7UDns2g) 

Portegies Zwart, S. F., Makino, J., McMillan, F., Hut, P. 2000b in prearation 

Spitzer, L. 1987, Dynamical evolution of globular clusters, Princeton, NJ, Princeton 
University Press, 1987, 191 p. 

Scalo, J. M., 1986, Fund, of Cosm. Phys., 11, 



This preprint was prepared with the AAS M5t;X macros v4.0. 



